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Abstract 

The medium modifications in the properties of QCD cascades are consid- 
ered. In particular, the changes in the intrajet rapidity distributions due to 
medium-induced decoherence, collisional losses of cascade gluons and those 
of final prehadrons are analyzed. 
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1 Introduction 



One of the most exciting experimental discoveries at RHIC are the dramatic 
changes in the patterns of energy flow in the events with large transverse 
energy. The simplest and most discussed related effect is the so-called jet 
quenching, the attenuation of hadrons with moderated and large transverse 
momenta, see the recent experimental review [T] and the theoretical reviews 
[21 m [31 El El [7] devoted to this problem. Generically the events with large 
transverse energies are jetty, i.e. are characterized by coUimated bunches 
of hadrons. A description of intrajet properties of high energy jets in usual 
QCD vacuum [S] belongs to the most successful applications of perturbative 
QCD. The detailed picture of parton cascade dynamics including fine details 
such as a destructive quantum interference leading to an angular ordering in 
sequential gluon decays was worked out and shown to successfully explain 
experimental data on intrajet spectra and multiplicity distributions. This 
explains a great interest in studying the modifications of intrajet properties of 
QCD jets produced in high energy nuclear collisions in which QCD cascades 
develop in high density environment. Experimental studies of this problem 
have been performed at RHIC [9j and are planned at LHC ^lOj . 

On the theoretical side, in order to study the medium-induced modifica- 
tions, one needs to work out the corresponding generalization of the formal- 
ism of the in- vacuum QCD. The presence of dense medium could influence 
the evolution of a cascade in two ways. First, the presence of the medium 
facilitates cascading by creating additional phase space for gluon emission. 
Second, there could be a non-radiative energy loss (gain) due to elastic in- 
teractions of the projectile with the modes of the medium. To account for 
interactions of the perturbative degrees of freedom in the evolving jet with the 
nonperturbative soft modes and/or medium one has to introduce an explicit 
mechanism of interaction between the partonic cascade and the environment. 

The most detailed studies of in-vacuum QCD cascades are performed 
with the help of the Monte-Carlo simulations. The two most popular ver- 
sions of these MC generators are PYTHIA [H] and HERWIG [12]. Both of 
them implement a key property of in-vacuum QCD jets, the angular order- 
ing, that in the modern versions is directly built in, see a detailed discus- 
sion in [16] and in the Appendix of the present paper. The first PYTHIA- 
based MC program including the effects of medium-induced radiative energy 
loss was PYQUEN [13]. Recently there appeared several MC generators 
treating the medium-induced effects in the QCD cascades, in particular the 
PYTHIA-based JEWEL [H] and Q-PYTHIA pjj and the HERWIG-based 
"Q-HERWIG" [H]. In the present paper we have used the object-oriented 
code realizing the 'global' 'constrained' mass-ordered version of PYTHIA im- 
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plemented in JEWEL [2]. All of the above-listed MC generators take into 
account the medium-induced radiative energy loss that is believed to be a 
major source of the energy loss; JEWEL, in addition, includes the effects of 
rescattering of cascade particles on those in the medium. The generic effect 
due to the medium-related interactions is the softening of intrajet particle 
spectra accompanied by the transverse broadening of jets. 

Let us note that for the case of non-radiative energy loss a correspond- 
ing modification of Altarelli-Parisi equations, constructed in analogy with a 
theory of ionization losses in electromagnetic cascades [T7], was considered 
in [iHl [19] where spectra and multiplicity distributions of QCD jets in a 
dissipative medium were studied. 

The focus of the present paper is on the analysis of some medium-induced 
modifications in the evolution of QCD cascades related to the fine structure 
of the spatiotemporal pattern of the in-medium QCD cascade, in particular 
on the substantial softening of the intrajet rapidity distribution. 

In the paragraph 12.11 we describe the basic features of the model of in- 
medium QCD cascade used in our simulations. 

In the paragraph 12. 21 we give the detailed description of the simulation of 
the mass-ordered QCD cascade. 

In the paragraph 12.31 we discuss the substantial softening of the intrajet 
rapidity distributions due to decoherence effects leading to the disappearance 
of the angular ordering. 

In the paragraph 12.41 we consider the effects of the collisional energy loss 
as related to the spatiotemporal pattern of the in-medium QCD cascade. 

We conclude the paper with the short summary of the results obtained 
and of the problems to be addressed in the near future. 

2 Decoherence and energy loss in in-medium 
QCD cascade 

In the present paper we will concentrate on discussing the properties of the 
timelike in-medium QCD cascades. The physical origin of the cascading 
process is the intense gluon radiation of a high energy parton created in a 
hard subprocess. Operationally one can think of this radiation cascade as of 
equipping the initial parton with energy Eq with some invariant mass Qq that 
is shaken off through decays into partons with lower invariant masses until 
reaching some non-perturbative scale Q/i/2 at which preconfined colorless 
clusters are formed. 

When considering the properties of the in-medium QCD cascades one has 
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to deal with several mechanisms modifying the in-vacuum intrajet charac- 
teristics. These are, in particular, the loss of quantum coherence and the 
resulting loss of angular ordering, collisional and radiative energy losses. In 
the present paper we shall focus on studying the impact of decoherence and 
collisional energy loss of partons and final prehadrons. The main motiva- 
tion is to analyze the dependence of observable quantities such as intrajet 
particle distributions, jet energy and multiphcity on the spatiotemporal pat- 
tern of the cascade, in particular on the distance from the origin of the 
cascade to the border of the hot dense fireball or, equivalently, the time in- 
terval during which the development of the cascade takes place inside the 
medium. To get a feeling for the corresponding physics, let us consider, for 
example, a gluon with the energy Eq = 100 GeV equipped with a typical 
initial invariant mass Qq = 10 GeV. This gluon will on the average decay 
at To ^ Eq/Qq = 1 GeV~^ ~ 0.2 fm, so that if the vertex in which the ini- 
tial hard gluon was generated was sufficiently far (several fermies) from the 
surface of hot fireball, a significant part of cascade vertices will be generated 
inside the medium. 

2.1 The model 

The key elements of the model developed in the paper are illustrated in Fig.[TJ 



The area to the left of the vertical dashed line is the region filled with 
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Figure 1: A sketch of the in-medium QCD cascade. 
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dense parton matter where modifications of cascade properties take place. 
To tfie riglit of tfiis line is the region of usual vacuum in which one has usual 
unmodified QCD cascade. The cascade is assumed to be initiated at the 
distance L from the medium-vacuum interface surface. 

The considered modifications of cascade properties inside the medium 
are related to several mechanisms: decoherence due to color randomization, 
collisional losses of cascading gluons and those of final prehadrons. 



2.1.1 Space-time pattern. 



To determine whether a given line is still inside the medium or has already 
left it, one needs to set up a clock counting time along the cascade. The only- 
possibility of associating a spatiotemporal pattern with the cascade is to use 
the lifetime r of a virtual parton which, for a parton with the energy E and 
virtuality that has been created in the decay of its parent parton with 
the virtuality QLr' reads 



1 



Ql 



par 



The lifetime of the initial parton is taken to be Tin = Eij^/Qf^. For example, 
the time Th at which the parton with the momentum qi in Fig. [T] produced 
the final prehadron is 



Ql 



l + z 



Ql 
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where z = Ei/Eq. For the parton under discussion it happened that Th < L 
so that both its formation at time Eq/Qq and its decay into final prehadron 
took place inside the medium. 



2.1.2 Decoherence due to color randomization 

For decays taking place inside the medium such as those of gluons with the 
momenta q and qi the angular ordering is absent. This means that, e.g., 
the magnitude of the angle 62 is not restricted by that of the angle 61. The 
corresponding effects will be discussed in the paragraph 12. 3[ 



2.1.3 Collisional losses of cascading gluons 

The gluons propagating in the medium (solid lines) experience, before con- 
verting to final prehadrons (black circles) collisional energy loss. The corre- 
sponding effects will be considered in the paragraph 12.41 
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2.1.4 Collisional losses of final prehadrons 

The final prehadrons formed inside the medium, such as, e.g., that formed by 
the gluon initially having the momentum qi, experience the medium- induced 
collisional loss and, if having insufficient energy, do not leave the medium. 
The effects related to the prehadron collisional energy loss are discussed in 
the paragraph 12. 4[ 

2.2 Mass-ordered QCD cascade 

In this paragraph we describe the procedure of generating a mass-ordered 
gluon cascade we employ. 

The cascade evolution goes through a sequence of decays q =^ qi + q2, 
where q = {E, is a four-momentum of the parent gluon and gi(2) are those 
of the daughter ones. We study the timelike evolution in which cascading 
reshuffies the initial jet virtuality in such a way that at each decay one has 

> 5'i + below. The cascading process at a given branch stops when 

a virtuality of the last parton reaches a threshold level Qi, i.e. no decay 
into two partons with minimal virtuality can occui0. The key property 
of the in-vacuum QCD cascades, effective angular ordering of subsequent 
decays due to color coherence, is not automatically taken into account in the 
mass-ordered evolution and has to be enforced explicitly by accepting only 
those generated decays for which the condition of angular ordering does hold. 

As the angular ordering is checked only a posteriori, the restrictions im- 
posed on initially generated decays are purely kinematical so that for the 
decay q qi + q2 the splitting variable z = Ei/E = l — E2/E should lie in 
the interval 

z4E,Q'\ QlQl) < z < hiE,Q'\ QlQl) (2) 

where 



z^{E,Q'\Ql,Ql)) = 




(3) 

Q± = V^Of i \/Q2- The restriction follows from the requirement of 
positivity of the square of transverse momentum generated in the decay. Let 
us also note that, as is immediately evident from (|3]), one has Qi + Q2 ^ Q- 
From equations (12|3|) it follows that the full specification of a decay re- 
quires, generally speaking, knowledge of invariant masses of the offspring 
partons. In our simulations we employ the so-called 'global' 'constrained' 

"'^By default the virtuality of such a last parton is set to 
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prescription that uses both the exact restriction fl2|5]) and its simphfied ver- 
sion in which Ql 2 are replaced by the minimal invariant mass Q\/ 4: 

z.{E,Q'\Ql)<z<z^{E,Q'\Ql) (4) 

where 



z^{E,Q'\Ql)^~z^ (e,Q'\^,^ 




h 

Q\ 

(5) 

In particular, the simplified restrictions (l4f5|) are used when calculating 
the Sudakov formfactor 



^(g2|E,QLx;QD = exp 



?2 

t max 



'4 r'""'"''d.°-'--<^-->''iiv.,/.w 

t Jz.{E,fi\Ql) 27r 



(6) 

The equation (E]) defines the probability to split into two gluons with virtu- 
alities Qhf^ for a gluon with the energy E having a virtuality restricted 
by condition Q\ < Q"^ < Q^ax- In ([6]) Ai'c = 3 is the number of colors and 

K{z) = — — — 2 + z{l — z) is the gluon splitting function. 

z[l — z) 

The explicit procedure we employ to generate an in-vacuum mass-ordered 
gluon cascade is as follows: 

1. First one draws, using (jHOnD, the scale at which the gluon under 
consideration branches into two new gluons. Let us note that at this 
step one fixes the lifetime of the gluon r = E{l/Q'^ — l/Qpar)- 

2. One draws the value of the splitting variable z determining the energies 
of the offspring gluons Ei = zE and Z2 = {1 — z)E. 

3. With the energies of the offspring gluons fixed, one draws their final 
invariant masses Q\ 2 and, therefore, fix their lifetimes 

^1,2 = £'(1,2) 



(1,2) 



The values of Q\ 2 are accepted if they do not violate ([3]). In other case 
the maximum of Qi 2 is redrawn until the condition ([3]) is satisfieci^. 



^Let us note that such a procedure tends to exclude particles that simultaneously have 
big energies and virtualities. Indeed, let us assume that we have drawn some Q E. 
This immediately leads to z w 1/2 and, therefore, makes it impossible to get Ei k, E with 
Qi K, El. At the same time if we drew Qi — Q — Qh/"^ and Q2 = Qh/'^^ we would get the 
case from ([5]): z « 1 — Qh/E. 
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4. Finally, one ensures angular ordering by drawing the splitting variables 
for the decays of the offspring gluons Zi^2 and accepting them if they 
satisfy the inequalities 



cor r esp ondingly. 

2.3 Decoherence of in-medium QCD cascades 

We have already mentioned that quantum coherence and resulting angu- 
lar ordering of gluon emissions plays a crucial role in the physics of QCD 
jets. The violent environment created in ultrarelativistic heavy ion colli- 
sions acts as a source of random energy-momentum and color with respect 
to in- vacuum processes. This kind of random impact tends to ruin the phase 
tuning lying at the heart of quantum interference phenomena. The situations 
in which random external influences destroy the interference-related effects 
are very common in solid state physics. A well-studied example is provided 
by studying the physics of weak localization in the presence of random ex- 
ternal impacts, see e.g. |22l [23]. Thus it is intuitively quite clear that the 
quantum coherence effects should be broken in the violent environment cre- 
ated in heavy ion collisions. This was first explicitly mentioned in [13], in 
which the origin of decoherence and ensuing disappearance of angular order- 
ing was related to rescattering of cascading particles on scattering centers in 
the medium. 

It is well known that in plasma physics the mean-field effects are typi- 
cally much stronger than those related to scattering. It therefore of interest 
to point out to another source of decoherence, the randomness of the color 
of mean field. A quantitative argument is provided by comparing the char- 
acteristic times of color and momentum diffusion in dense strongly interact- 
ing matter, tc and tp correspondingly, that were calculated for quark-gluon 
plasma in pT| : 



where is an electric screening mass and rriM - 9. (non-perturbative) mag- 
netic screening one. Obviously the pace of color randomization is much 




tp « [4alT\n{l/as)] 
tc ^ [3asTln{mE/mM)] 



-1 



(7) 
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faster than that of momentum evolution and, therefore, coUisional energy 
loss, tc tp- This situation is similar to decoherence related to rotation 
of spin by random external magnetic filed [23]. Therefore, when describing 
a propagation of a test high energy colored particle through plasma, one 
can assume that at time scales of order of tc all coherent effects related to 
precise color matching are destroyed. The origin of angular ordering in the 
in-vacuum QCD cascades is exactly in such kind of precise color matching, 
therefore at the timescales r > tc the color coherence is broken. Of course, 
the estimates for the timescales in ([7j) refer to equilibrated QCD matter con- 
sidered in the hard loop approximation and are thus not too realistic. In the 
context of the present paper it is only important that ([7]) points out to color 
randomization as the fastest source of decoherence in the hot dense fireball. 

In what follows we shall use a simplest assumption that angular ordering 
is always broken for gluon decays in QCD cascades that take place inside the 
hot dense medium. In this paragraph we consider the situation in which such 
breaking of angular coherence is the only medium-related effect in the gluon 
cascade. It is known for a long time that in the absence of angular ordering 
the properties of QCD cascades are very different from those when angular 
ordering is imposed [25]. In particular, the jet multiplicities decrease, but the 
strongest effect is related to the changes in intrajet energy distribution which 
experiences dramatic softening [25]. As in our considerations the effects of 
breaking of color ordering are present only for decays taking place inside the 
medium, the cumulative effect on energy distribution of final prehadrons, 
the strength of the above-mentioned softening should depend on the relative 
proportion of vertices generated inside the medium or, in other terms, on the 
distance L from the cascade origin to the border of the fireball. 

In Fig. [2] we plot the distributions of final prehadrons in the rapidity 
y = ln(£'o/-E') for gluon jets with the initial energy Eq = 100 GeV and L = 
0.5, 1 and 5 fm. 

From the distributions in Fig.[2]we see that the effect of color decoherence 
is indeed very strong. Already at L = 0.5 fm the distribution is significantly 
different from the vacuum one and at larger L practically saturates. Thus, 
even in the absence of energy loss in terms of energy-momentum, the intrajet 
properties are significantly affected by the medium through color decoher- 
ence. 

The changes in distribution over the rapidity \n{E^et/E) shown in Fig. |2] 
reflect important changes in the spatiotemporal pattern of the jet. Indeed, 
from the inside-outside nature of particle production in which soft particles 
are produced faster than hard ones, cf. Eq. ([T]) one can conclude that the 
substantial softening of the rapidity distribution shown in Fig. |2] reflects an 
increased yield of short formation times of prehadrons. Thus it is of interest 
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log(Ejet/E) 

Figure 2: Rapidity distribution P{y) of final prehadrons: 1) L = Ofm, full 
angular ordering, red, solid; 2) L = 0.5 fm, partial angular ordering, green, 
dashed; 3) L = 5 fm, partial angular ordering, blue, dotted. 

to consider the yield of prehadrons formed inside the medium of given length. 
This yield is plotted in Fig. [3] for jets with and without angular ordering. 

2.4 Collisional losses in the in- medium QCD cascades 

In this paragraph we shall make a crude estimate of the efects of the part 
of medium-induced energy loss which is proportional to the distance / that 
a parton or prehadron travel inside the medium. For simplicity, in analogy 
with QED cascades in matter, we shall term these losses as collisional losses. 
Radiative medium-induced losses are, at least at small times, proportional to 
and require a separate treatment. Usually one considers the radiative losses 
as dominant, so our restriction to the collisional ones provides an estimate 
of the minimal medium-induced energy loss. 

In analogy with QED showers [2SIEZ], the partonic intrajet energy losses 
AEp are described by the continuous distribution V{AEp, l\ fi^E, <^ae)- The 
parameters {fiAE, o'ae) describe the mean energy loss and the corresponding 
standard deviation at the distance of 1 fm. In our simulations we use a simple 
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Figure 3: Relative yield of prehadrons formed inside the medium: 1) global 
angular ordering, red, solid; 2) partial angular ordering, green, dashed; 3) 
partial angular ordering and collisional losses, blue, dotted. 



Gaussian parametrization of V{AE, l\ /iae, o'ae)'- 

V{AEp,l\ fiAE,aAE) = uS{AEp)+-=i -exp Q{AE,) 

(8) 

where 

1 /"°° 1 
u = -j= / dx exp{-x^ /2) = -(1 - erf(V2/2)) ^ 0.159 
v2n Ji 2 

as we use only the case /iae = ctae- The numerical values we chose in most 
of the simulations were /iae = 1 GeV/fm and cr^^ = 1 GeV^/fm. 

The losses AEh experienced by prehadrons formed inside the medium on 
their way out of the hot fireball were simply taken to be proportional to the 
distance they cover inside the medium, AE^ = 1 GeV/fm ■ /. 

Taking into account the softness of the majority of produced prehadrons 
that takes place even without taking into account the energy losses, just 
because of decoherence, an interesting question to address is how many pre- 
hadrons have enough energy to leave the medium when the energy losses 
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Figure 4: Relative yield of particles stopped inside the medium: 1) parton 
losses, red, solid; 2) parton and hadron losses, green, dashed. 



are switched on. The simplest crude criterion of stopping one can use is 
Ecrit = Qh/'^ SO that the momentum of the corresponding parton or pre- 
hadron is zero. The corresponding yields are shown in Fig. |H 

We see that at large enough distances to the border of the fireball the 
effect is very significant: at L = 5 fm more than half of the final prehadrons 
are stopped inside the fireball. The effect of stopping is also clearly seen in 
Fig. [5] in which we plot the average multiplicity of final prehadrons inside the 



cone 6 



jet 



0.7. 



For completeness in Fig. [H] we show the dependence of average energy 



within the cone of 6-, 



jet 



0.7 on L. 



The additional softening incurred by energy losses does considerably change 
the shape of the rapidity distributions of final prehadrons. In Fig. [7] we plot 
the corresponding distributions at L = 5 fm and in Fig. [8] we illustrate the 
dependence of rapidity distribution taking into account all sources of energy 
loss on L. 
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Figure 5: Mean jet multiplicity as a function of medium size: 1) no losses, 
red, solid; 2) intracascade parton losses, green, dashed; 3) intracascade and 
final state losses, blue, dotted. 
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Figure 6: Mean jet energy as a function of medium size: 1) no losses, red, 
solid,; 2) intracascade parton losses, green, dashed; intracascade and final 
state losses, blue, dotted. 
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Figure 7: Rapidity distribution P{y) of final prehadrons, cascade with partial 
angular ordering: 1) no losses, red, solid; 2) intracascade parton losses, green, 
dashed; 3) intracascade and final state losses, blue, dashed. 
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Figure 8: Rapidity distributions at different medium length L 
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3 Conclusions 



The main goal of the present study was to analyze some medium-induced 
changes in the the evolution of the in-medium QCD cascade. 

In the focus of our attention was a detailed analysis of the medium- 
induced color decoherence and the ensuing disruption of angular ordering for 
the decays taking place inside the medium. We have shown that this effect 
leads to the substantial softening of the intrajet rapidity spectrum. 

Another quantitatively important effect considered in the paper is that of 
the coUisional losses of cascade gluons and of those final prehadrons that hap- 
pened to be formed inside the medium. We have shown that, at variance with 
the widespread opinion that the dominant majority of the final hadrons form 
outside the hot and dense fireball, the yield of prehadrons formed inside the 
fireball is quite significant. Combined with the effect of partonic collisional 
loss, this leads to the substantial deformation of the rapidity spectra. 

To make the analysis presented in this paper more complete one has, of 
course, take into account the medium-induced radiative losses. Work on this 
problem is currently in progress. 
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Appendix 

Let us discuss the approximations which made when "angular orderung" 
conditions is introduced. Recall that the "angular" vari bid which is used in 
Herwig and similar is written as 

EbEc Ei,Ec 2 EbEc 

where we have used 

fa = (Pb + Pcf=pt+fc+ '^{Pb, Pc) (A.2) 

The initial upper limit for ^ usually supposed as 1. This means that we 
consider the openning angles 9a to be less then tx/2 only. In this area the 

■^This is exactly the same variable for which the ordering was originally derived in [2 8) 
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conditions of angular ordering (9) and sin 6' (or pt or sin^ 6') are equivalent. 
In terms of momenta we have for sin^ 9 



sm 



1 t/a = i- ^ ^ = i T ZTT^^ 

where we again used f lA.2p . 

Pb 



(A.3) 





Figure 9: Angular ordering 

As for considered angles pi — pi — pI > Owe can rewrite the ordering 
condition as following (see Figure for parton indexes) 



-Q -0 -Q -fi -fi -fl 

Pc-Ps-Pd ^ Pa-Pb- Pc 



\Ps\\Pd 



\Pb\\Pc 



and for ^ variable 



-Q -0 -Q -0. -fl -fl 

Pc-Ps- Pd ^Pa-Pb-Pc 



(A.4) 



(A.5) 



EgEd EbEc 
Comparing (1A.4P and (lA.SI) we see that ^-ordering coinside with true angular 



ordering if we suppose E ^ \p\ in denominators of (]A.4p . Which condition 
is stronger depends on the four ratios Q/E and can be different for different 
cases. 

From other side 



^ ^-El-El-El-Ql + Ql + Ql = 2E,E^-Ql + Ql + Ql{k.Q) 



Pa-Pb-Pc 



where we have used energy conservation Ea = Eh + E^. And we can rewrite 
condition (lA.Sp as following 



Ql-Qi-Ql^Qi-Qi-Qi 



Eg Ed 



EbEc 



(A.7) 
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E, 

or taking into account that EsE^ = Zc{l — Zc)E'^ and Eh = ZaEa = 



I - z„ 



Ql-Ql-Ql ^ {Ql-Ql-Ql){i-Za) 

If we now suppose that > + Ql and > + we get the 
conditions which is used in "Virtuahty ordered" Pythia as angular ordering 
condition 

^ « (A.9) 
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